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Abstract. Particle acceleration in massive star forming regions can proceed via a large 
variety of possible emission scenarios, including high-energy gamma-ray production in 
the colliding wind zone of the massive Wolf-Rayet binary (here WR 20a and WR 121a), 
collective wind scenarios, diffusive shock acceleration at the boundaries of wind-blown 
bubbles in the stellar cluster, and outbreak phenomena from hot stellar winds into the in- 
terstellar medium. In view of the recent Fermi-LAT detection of HESS J1023-575 (in the 
vicinity of Westerlund 2), we examine another very high energy (VHE) gamma-ray source, 
HESS J1848-0145 (in the vicinity of W43), possibly associated with a massive star cluster. 
Considering multi-wavelength data, in particular TeV gamma-rays, we examine the avail- 
able evidence that the gamma-ray emission coincident with Westerlund 2 and W43 could 
originate in particles accelerated by the above-mentioned mechanisms in massive star clus- 
ters. 
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1. Introduction 

Young stellar clusters and star forming re- 
gions have long been proposed as gamma-ray 
sources. This relates to the existence of late 
stages of stellar evolution, which are gener- 
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ally considered favorable for particle acceler- 
ation. Different scenarios have been suggested 
for the gamma-ray production in these environ- 
ments, considering the contribution of hot and 
massive Wolf-Rayet (WR) and OB stars and 
their winds in single, binary, or collective pro- 
cesses, pulsars and their synchrotron nebulae, 
as well as supernova explosions and their ex- 
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panding remnants (Reimer et al. 2006; Romero 
2010). Their detection and identification us- 
ing gamma-ray data have already a rather long 
story, first attempts being made with COS- 
B (Montmerle 1979) and later with CGRO- 
EGRET (Romero et al. 1999). 
Several years after EGRET, the Large Area 
Telescope (LAT), aboard Fermi, offers the op- 
portunity to study faint and extended gamma- 
ray sources such as star forming regions, mas- 
sive star clusters and pulsar wind nebulae 
(PWNe). The LAT is an electron-positron pair 
conversion telescope, sensitive to y-rays with 
energies between 30 MeV and 300 GeV, with 
improved performance (a large effective area, 
a broad field of view, and a very good an- 
gular resolution) compared to its predeces- 
sor (Atwood et al. 2009) . 
The morphological and spectral analyses pre- 
sented below were performed using two dif- 
ferent methods to validate our results: a 
maximum-likelihood method (Mattox et al. 
1996) implemented in the Fermi SSC sci- 
ence tools as the "gtlike" code and "point- 
like", a binned likelihood technique that op- 
erates by convolving the intrinsic LAT point 
spread function (PSF) with a spatial model of 
the candidate source to create a "Pseudo-PSF" 
of the source appearance (Kerr 2010). 

2. Westerlund 2 

The prominent giant H II region RCW 49 and 
its ionizing cluster Westerlund 2 are located 
towards the outer edge of the Carina arm of 
our Milky Way at a distance uncertain in the 
range of values between 2.2 kpc (Brand & Blitz 
1993) and 8 kpc (Rauw et al. 2007). RCW 
49 is a luminous massive star formation re- 
gion, which hosts over two dozen of massive 
stars as well as two remarkable WR stars (WR 
20a and WR 20b), and has been studied ex- 
tensively over the whole electromagnetic spec- 
trum. This interest was further boosted by the 
HESS discovery of the bright and extended 
VHE gamma-ray source HESS J1023-575 in 
the vicinity of Westerlund 2 (Aharonian et al. 
2007). The extent of the observed gamma-ray 
emission spatially overlaps with Westerlund 2 
and the WR stars WR 20a and WR 20b. Even 



though the source extension does not support 
a colliding wind scenario or an emission pro- 
duced by the Westerlund 2 cluster itself, other 
scenarii include diffusive shock acceleration 
in the wind-blown bubble itself and super- 
sonic winds breaking out into the interstellar 
medium. 

Recently, Fermi-LPtf discovered the very 
young (characteristic age of 4.6 kyr) and en- 
ergetic (spin-down power of 1 . 1 x 10 37 erg s" 1 ) 
pulsar PSR J 1023 -5746, coincident with the 
TeV source HESS J1023-575 (Saz Parkinson 
et al. 2010). This re-opened the ques- 
tion of the identification of the gamma-ray 
source and led to re-observations and anal- 
yses by HESS (Abramowski et al. 2010) 
and Fermi (Ackermann et al. 2011). A sig- 
nificant signal in the off-pulse window of 
PSR J 1023 -5746 was obtained using Fermi 
data, with an emission above 10 GeV signif- 
icant at more than 3<x level. In addition, the 
GeV flux connects well with the one derived 
at TeV energies (as seen in Figure 1) with 
a luminosity corresponding to less than 1.5% 
of the pulsar spin-down (assuming a distance 
of 2.4 kpc). This supports highly the scenario 
in which the gamma-ray source detected by 
Fermi and HESS would be the PWN pow- 
ered by PSR J1023-5746. On the other hand, 
the molecular content of the region provides 
sufficient target material to explain the emis- 
sion through hadronic interaction of cosmic- 
rays accelerated by a range of feasible mech- 
anisms in the open cluster interacting with 
molecular clouds. In addition, the lack of X-ray 
emission from the immediate vicinity of PSR 
J1023-5746 (Fujita et al. 2009) is perplexing 
given its extremely high spin-down luminos- 
ity. More data especially in X-rays and TeV are 
therefore needed to elucidate the real nature of 
HESS J1023-575. 

3. W43 

The W43 star forming region is located along 
the Scutum-Crux spiral arm tangent at a dis- 
tance of ~ 6 kpc (Nguyen Luong et al. 2011). 
The core of W43 harbors a well-known giant 
H II region powered by a particularly lumi- 
nous cluster of WR and OB stars, including 
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Fig. 1. Spectral energy distributions of the off-pulse 
emission of PSR J1023-5746 (Ackermann et al. 
2011). The LAT spectral points (red) are obtained 
using the maximum likelihood method "gtlike" in 
7 logarithmically-spaced energy bins. A 95 % C.L. 
upper limit is computed when the statistical signif- 
icance is lower than 3 cr. The blue points represent 
the H.E.S.S. spectral points (Aharonian et al. 2007). 
The Suzaku upper limit is shown with a green arrow 
(Fujita et al. 2009). The black line denotes the total 
synchrotron and Compton emission from the nebula. 
Thin curves indicate the Compton components from 
scattering on the CMB (long-dashed), IR (medium- 
dashed), and stellar (dotted) photons. 



WR 121a (a WN7 subtype star), which is char- 
acterized by extreme mass loss rates. Motte 
et al. (2003) established it as a dense region 
equivalent to a mini-starburst region since it is 
undergoing a remarkably efficient episode of 
high-mass star formation. Recent far-infrared 
to submillimeter data from the Herschel Space 
Observatory revealed a complex structure of 
chimneys and filaments, and confirmed its ef- 
ficiency in forming massive stars (Bally et al. 
2010). A detailed 3D (space-space-velocity) 
analysis of the molecular and atomic cloud 
tracers through the region has been carried, re- 
vealing that W43 is particularly massive and 
contains a total mass of 7.1 x 10 6 Mq of 
molecular gas at a distance of 6 kpc (Nguyen 
Luong et al. 2011). 

At very-high energies, HESS reported a signif- 



icant excess (9<x pre-trials) with an extended 
and complex morphology (Chaves et al. 2008). 
The TeV source HESS J1848-018 is in the 
direction of, but slightly offset from, the star 
forming region W43 and there are no ener- 
getic pulsars or supernova remnants within 
0.5° which seems to favour an association with 
the star forming region W43. 
Trying to unveil the nature of the TeV source, 
we analyzed 3 1 months of Fermi data collected 
starting August 4, 2008, and extending un- 
til March, 14, 2011. Only gamma-rays in the 
Pass 7 Source class events were selected, and 
from this sample we excluded those coming 
from a zenith angle larger than 100° because 
of the possible contamination from secondary 
gamma-rays from the Earth atmosphere (Abdo 
et al., 2009d). Figure 2 shows a TS map of 
the region in the energy range above 2 GeV: a 
significant emission coincident with the HESS 
source and slightly offset from WR 121a is 
clearly visible. Since the TeV source is sig- 
nificantly extended, we tried to determine the 
extension of the Fermi source using pointlike: 
a 3.7cr improvement is obtained by using a 
Gaussian width of 0.3° in comparison to a sim- 
ple point source. However, the statistics is not 
large enough to discriminate between one ex- 
tended source and several point sources. The 
spectral analysis performed above 100 MeV 
assuming the best Gaussian model shows that 
the Fermi source is well described by a log- 
parabola model with spectral parameters a 
2.07 + 0.07 stat , /3 = 0.22 + 0.05 stat , E b = 1.0 
GeV, and an integral flux above 100 MeV of 
(1.77 + 0.26 stat ) x 10~ 7 cm" 2 s' 1 : 



dN _ IE 
dE I Et, 



-(a+/3\og(E/E b )) 



(1) 



The SED presented in Figure 3 is very similar 
to that obtained on most pulsars detected 
by Fermi-LAiT which could indicate that we 
are contaminated by the emission coming 
from one energetic radio-faint pulsar (since 
there are no radio-loud pulsars in this region), 
as it is the case for Westerlund 2. A blind 
search is on-going: a pulsed detection in 
this extended region is a real challenge (the 
extent of our source is 0.3°): indeed, when 
the pulsar is off by only 1 arcmin it is already 



Lemoine-Goumard M. et al.: Fermi-LAT observations of W43 and Westerlund 2 



285 



iLllS50.7-001*F'jLJlS19.3-005S 



2FQLJ1849 9^^H ' J 



2FGLJ1847.2-0236 

X 



32.500 32.000 31.500 31.000 30.500 30.000 



Fig. 2. Test Statistic (TS) map of the region sur- 
rounding W43 with side-length 3°, above 2 GeV. 
Each pixel of this image contains the TS value 
for the assumption of a point-source located at the 
pixel position, Galactic diffuse emission and nearby 
sources being included in the background model. 
The position and extent of the HESS source is repre- 
sented with a white circle. The position of WR 121a 
is indicated with a green cross. Sources extracted 
from the 2FGL catalog and taken into account in the 
analysis are marked with a blue cross. 



hardly detectable ! As for Westerlund 2, 
multi-wavelength data would greatly help to 
reveal the nature of the gamma-ray source, 
especially in X-rays and in TeV. 

These two examples indicate that stel- 
lar clusters and massive star formation re- 
gions might be gamma-ray sources. However, 
this identification requires first a discrimina- 
tion with the Pulsar/PWN (GeV/TeV) scenario 
which is a real challenge in these crowded re- 
gions of the Galactic Plane. 

Acknowledgements. The Fermi LAT Collaboration 
acknowledges support from a number of agencies 
and institutes for both development and the oper- 
ation of the LAT as well as scientific data analy- 
sis. These include NASA and DOE in the United 
States, CEA/Irfu and IN2P3/CNRS in France, ASI 
and INFN in Italy, MEXT, KEK, and JAXA in 
Japan, and the K. A. Wallenberg Foundation, the 
Swedish Research Council and the National Space 





j 


PRELIMINARY 






Statistical errors or! 


10" 


i Fe 










■c - 






■ 


L' L L1111IL 





n: 1 



1<£ 

Energy [MeV] 



Fig. 3. Spectral energy distribution of the source 
coincident with HESS J1848— 018 obtained using 
pointlike. Only statistical errors are shown. The red 
solid line presents the best fit obtained assuming 
a log-parabola. The HESS spectrum (Chaves et al. 
2008) is indicated with the green solid line. 
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